Draft version February 5, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 08/13/06 



THE ROTATING MOLECULAR CORE AND PRECESSING OUTFLOW OF THE YOUNG STELLAR OBJECT 

BARNARD IC 

Brenda C. Matthews 

Herzberg Institute of Astrophysics, National Research Council of Canada, 5071 West Saanich Road, Victoria, BC, V9E 2E7, Canada 



MiCHIEL R. HOGERHEIJDE 
Leiden Observatory, P.O. Box 9513, 2300 RA, Leiden, The Netherlands 



JeS K. J0RGENSEN 

^ 1 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street MS42, Cambridge, MA, 02138 

! AND 

Edwin A. Bergin 

^ ^, Department of Astronomy, University of Michigan, 825 Dennison Building, 500 Church Street, Ann Arbor, MI 48109 

^ , Draft version February 5, 2008 

m ■ 

; ABSTRACT 

■ We investigate the structure of the core surrounding the recently identified deeply embedded young 
stellar object Barnard Ic which has an unusual polarization pattern as traced in submillimeter dust 

, emission. Barnard Ic lies within the Perseus molecular cloud at a distance of 250 pc. It is a deeply 

^ ■ embedded core of 2.4 Mq (Kirk et al.) and a luminosity of 4 ± 2Lq. Observations (and resolutions) 

O ' of ^^CO J = 1 - (9.2" X 5.9"), i^co J = 1 - and C^^O J = 1 - (14.3" x 6.7"), HCO+ 

0> ! J = 1 - (7.6" X 5.8"), and N2H+ J = 1 - (5.9" x 4.6") were obtained with the Berkeley-Ilhnois- 

10 ' Maryland Association array, together with the continuum at 3.3 mm (6.4" x 4.9") and 2.7 mm (9.5" 

0^ ' X 6.3"). The field of view of the BIMA array antennas at 3 mm is 2.1'. Single-dish measurements of 

, N2H+ J = 1 — and HCO"*" J = 1 — with FCRAO reveal the larger scale emission in these lines 

■ with resolutions of 57.5" and 60.5", respectively. The ^^CO and HCO+ emission traces the outfiow 
I extending over the full field of view, which coincides in detail with the S-shaped jet recently found in 

r~| . Spitzer IRAC imaging. The N2H+ emission, which anticorrelates spatially with the C^^O J = 1 — 

pi i' emission, originates from a rotating envelope with effective radius ^ 2400 AU and mass 2.1 — 2.9 Mq, 

as derived from the 3.3 mm continuum emission. N2H+ emission is absent from a 600 AU diameter 
region around the young star, offset from the continuum peak. The remaining N2H+ emission may lie 
in a coherent torus of dense material. With its outflow and rotating envelope, Blc closely resembles 
^ I the previously studied object L483-mm, and we conclude that it is a protostar in an early stage of 

• • , evolution, i.e., Class or in transition between Class and Class I. We hypothesize that heating by 

the outflow and star has desorbed CO from grains which has destroyed N2H+ in the inner region and 
surmise that the presence of grains without ice mantles in this warm inner region can explain the 
unusual polarization signature from Blc. 
' Subject headings: ISM: clouds — ISM: molecules — ISM: individual (Barnard 1) — stars: formation 

— radio lines: ISM 
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1. INTRODUCTION 

Class sources represent the youngest phase of 
low-mass star formation. They are characterized 
by higher infall rates than more evolved sources 
based on outflow activity (Bontcmps ct al. 199_^ 
iWhitworth fc Ward-Thompsonll2001l|l . an absence of op- 
tical emission, and a high ratio of submillimeter to bolo- 
metric luminosity. Observations of the internal struc- 
ture of Class objects rely on interferometers because 
the sources remain deeply embedded within their par- 
ent molecular clouds. To understand the kinematics 
of the collapse process requires observation of multi- 
ple molecular lines because the chemistry within these 
dense, cold cores is complex. Carbon-bearing species 
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are observed to be strongly depleted within the core in- 
teriors due to freeze-out onto grains. Nitrogen-bearing 
species were thought to deplete much more slowly than 
carbon-bearing species because of the low binding en- 
ergy of the N2 molecule, requiring long periods of time or 
very high densities to show depletion l)Bergin fc Langeil 
^97; Bergin et al. 2002). However, recent laboratory 
results indicate that CO binding energies are compa- 
rable to those of N2. These studies also indicate that 
the degree of mixing of N2 and CO ices impacts the 
desorption timescales fObcr g et al.ll2005t iBisschop et alJ 
i200 6'). Nonetheless, it is clear that N2H"'" (and its daugh- 
ter product NH3) show a strong rise in abundance as 
its main destroyer, ^^CO, depletes from the gas, mak- 
ing N 2H+ an effective tra cer of internal core kinematics 
Ce.g.. lAikawa et aLll2005j) . When ^^CO is present how- 
ever, reactions between C+ or ^^CO lead to the eventual 
destruction of the two standard nitrogen-bearing tracers: 
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NH3 and N2H+. 

Recent observations of several Class sources have 
shown an absence of N2H+ emission at the source po- 
sition on sc ales within several hundred AU of the pro - 
tostar fi.e.. I.Torgensen. Schoier fc van Dishoeck ll2004bD . 
Typically, N2H+ emission is depressed at the core center 
and has been interpreted as destruction of N2H+ due to 
evaporation of CO and its isotopes from dust grains (e.g., 
L483-mm (hereafter L483), Joracnscn 2004) or depletion 
of N2H+ at high densities (e.g., IRA M 04191+1 522 (here- 
after IRAM 04191), Bellochc, A ndre. Despois fc Bhnde^ 
12002(1 . An absence of N2H"'" cannot be purely an indica- 
tor of age; L483 is th ought to be in transition between 
Class and Class I l|Tafalla et al.l bOOO) while IRAM 
04191 is estimated to be among the youngest Class 
sources known CAndre. Mottc & Bacmann 19991. Inter- 
estingly, in both cases, the N2H+ emission morphology 
is double peaked around the source center with a signa- 
ture indicating rotation, and an anticorrelation is noted 
between N2H'^ emission and the emission from tracers 
dominated by the outflow (^^CO and HCO+). However, 
in L483, C^^O is centrally peaked on the source indicat- 
ing that ^^CO has been evaporated from dust grains at 
the core center. The combination of N2H''" and C^*0 
morphology could thus discriminate between whether an 
absence of N2H+ is due to depletion onto dust grains or 
the destruction of N2H+ in the presence of ^^CO. 

In this paper, we present observations of the inter- 
nal structure of the protostellar core Barnard Ic (Blc) 
from the Berkeley-Illinois-Maryland Association (BIMA) 
array and the Five Colleges Astronomical Radio Ob- 
servatory (FCRAO) 14 m telescope. Blc was discov- 
ere d during 850 Atm polarimet ry mapping in Barnard 1 
bv IMatthews fc WilsonI pOOl . Recent IRAC imaging 
from Spitzer reveals that this source is highly reddened 
and deeply embedded in the Bl clo ud with an exten- 
sive, highly colhmated outflow (Jo rgensen et alJl2006j) . 
The presence of central cavities within young protostellar 
cores is of particular relevance to Blc because a heated 
central region could explain why Blc has a unique sig- 
nature in polarized emission. Its polarization pattern 
suggests th at the polarized intensity rises to the center 
of the core ((Matthews fc Wilsonll2002D . rather than flat- 
tening out as seen in other cores. When compared to the 
total intensity, a flat distribution in polarized intensity 
produces a declining ratio toward the peaks of cores. All 
other low-mass starless and star-forming cores observed 
in polarized dust emission have a so-called "polarization 
hole" at high intensities, thought to arise from changes in 
magn etic field geometry or dust grain physics ijMatthewsl 
I2005|) . One of the favored explanations for the polariza- 
tion holes is that the grains at core centers are ineffec- 
tive polarizers (due to changing grain physics). Heat- 
ing removes the outer grain mantles and could increase 
the polarization efficiency of the grains within cavities 
(Whittet et al. 2001). 

Barnard 1 is part of the Perseus molecular cloud 
complex (pachiller fc Cernicharol [1986(1 which is one of 
the closest star-forming regions to the Sun. Its dis- 
tance is the subject of some debate, with estimates 
ranging from 200 pc (based on extinction studies, 

ICernisI IT990I) to 330-350 pc (based on the Pe r 0B2 
association, Ide Zeeuw. Hoogerwerf fc de Bruiin3 Il999t 



iBorgman fc Blaauwl[T96l iHerbig fc Jonejll983(l . It has 

been suggested that the compl ex may be comprised o f 
two clouds at 200 and 300 pc l|Cernicharo et al.lll985l). 
that there is a distance gradient along its length l|Sargentl 
I1979|) . or that it is comprised of independent clouds at 
varying distances [Ridge et al. 200 6). To facilitate com- 
parison with recent works on the Per seus cloud (i.e., 
lEnoch et alJ2OO6tlJ0rgensen et al.l2006j) . we adopt a dis- 
tance of 250±50 pc to the Barnard 1 cloud, as determined 
from recent estimates of extinction (|Cernis fc Straizv^ 
12003 ) and measure ments of parallax in members of IC 
348 l)Belikov et al.l l20021 . 

This paper presents high resolution intcrferometric 
data from several molecular species and continuum emis- 
sion at 2.7 and 3.3 mm. The main objective is to deter- 
mine whether Blc exhibits a central cavity which could 
help explain its unique polarization properties. The ob- 
servations and data reduction techniques are described 
in § 121 We present the continuum results and derive the 
core mass in § |2| The molecular line data are presented 
in § ^ We discuss these data in § |5l Our findings are 
summarized in § 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. BIMA Interferometric Data 

Observations were made over the period of 2002 Oc- 
tober to 2003 April using the Berkeley-Illinois-Maryland 
Association (BIMA) interferometer ( Welch ct al. 199(|) 
in Hat Creek, CA. We utilized a single pointing toward 
the position aj2ooo = 03''33"n7;8, 6j2ooo = -l-31°09'33'.'0 
(J2000). Two configurations of the ten 6.1 m antennas 
were used to observe the lines of N2H+ J = 1 — and 
HCO+ J = 1—0. The C-array and B-array had projected 
baselines between 2-33 kA and 3-74 kA, respectively. The 
12C0 J = 1 - line and its isotopomers ^^CO J = 1 - 
and C^^0J = 1 — were observed only in the C-array 
configuration. Tabled contains the sensitivities achieved 
per track. 

The N2H+ J = 1 — line was observed utilizing the 
digital correlator to record the line in bands of 6.25 and 
12.5 MHz width with 256 channels each, giving resolu- 
tions of 0.079 and 0.157 km s~^. The larger bandwidth 
window permitted detections over the range of the seven 
hyperfine components of N2H+ J = 1 — 0. All 10 anten- 
nas were available for both 93 GHz tracks. 

The HCO+ J = 1 — line was observed in bandwidths 
of 6.25 and 12.5 MHz with 256 channels each, resulting 
in resolutions of 0.082 and 0.164 km s~^, respectively. 
A window in the upper side band was senstive to SiO 
J = 2 - 1 in a 12.5 MHz window with 0.085 km s^^ 
resolution. Three tracks were obtained in the B array 
and two in the C array. All 10 antennas were available 
for three tracks, with one antenna offline for a B array 
track and two antennas missing from one of the C array 
tracks. 

The ^^CO J = 1 — line was observed in a band 12.5 
MHz wide across 256 channels. We also observed it in 
100 MHz windows to detect high velocity CO gas. The 
velocity resolutions were 0.127 km s~^ and 8.127 km s"'^, 
respectively. No CO emission was detected at velocities 
exceeding 20 km s^^ from the rest velocity of the source. 
As for the N2H+ J = 1 — datasct, the phase calibrators 
were 3c84 and 0237-1-288. Only 9 of 10 antennas were 
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TABLE 1 

Observational Summary of Molecular Species 



Tr3,iisition 


FrGC(UGiicy 


Date 




Resolution 


Sensitivity 








Configiircition 


[Kin h J 


[Jy beam 


N2H+ J = 1 - 


93.17378 


19 November 2002 


c 


0.079 


0.4 










0.157 


0.28 






27 February 2003 


B 


0.079 


0.32 










0.157 


0.21 


HCO+ J = 1 - 


89.18852 


24 October 2002 


C 


0.082 


0.4 










0.164 


0.27 






25 January 2003 


B 


0.082 


1.3 










0.164 


0.9 






11 March 2003 


B 


0.082 


1.1 










0.164 


0.75 






16 March 2003 


B 


0.082 


0.6 










0.164 


0.4 






07 April 2003 


C 


0.082 


0.38 



0.164 0.26 



SiO J = 


2 


- 1 


85.64046 


24 October 2002 


C 


0.085 


0.38 










25 January 2003 


B 




1.4 










11 March 2003 


B 




1.0 










16 March 2003 


B 




0.5 










07 April 2003 


C 




0.37 


12CO J = 


1 


- 


115.21720 


04 December 2002 


C 


0.127 


1.4 


"CO J = 


1 


- 


110.20135 


20 November 2002 


c 


0.066 


0.65 














0.133 


0.45 


C180 J = 


1 


- 


109.78217 


20 November 2002 


c 


0.067 


0.67 














0.133 


0.48 



''Not detected. 

useable for this track due to phase incoherence on an- 
tenna 9. We were sensitive to the CN J — 1 ~ line in 
the USB during this observation, but none was detected. 

The i^co J = 1 - and C^^O J = 1-0 hues were 
observed with the digital correlator configured to record 
the lines in the upper side band with 12.5 MHz band- 
width over 256 channels. One antenna was offline during 
this track. Data from a second antenna was flagged due 
to the same phase problems. Spectral line data from a 
third antenna was removed due to noisy phases and am- 
plitudes. 

Phase and amplitude variations were calibrated by ob- 
serving the nearby quasars 3c84 and 0237-1-288 (when 
3c84 reached elevations exceeding 85°) approximately ev- 
ery 30 minutes. The adopted fluxes of these quasars 
was epoch dependent and measured against observa- 
tions of the planet Uranus when possible. The cali- 
bration was performed using the MIRIAD (Multichan- 
nel Image Reconstruction, Ima ge Analysis and Display; 
ISault. Teuben fc Wrightl JlQOl 'l task MSELFCAL. Ab- 
solute flux calibration was done using Uranus when ob- 
served or by derived fluxes of the gain calibrators during 
the same epoch as our observations (utilizing the cata- 
logue of fluxes at "plot_swflux" on the BIMA website^). 
Based on the uncertainty in flux of the calibrator and the 
relative variations in the flux of the quasar, we estimate 
our overall flux calibration is accurate to the 30% level. 

Subsequent processing of the data, including the com- 
bination of data from different conflgurations, were done 
with MIRIAD. Images were produced using MIRIAD 's 

^ http://bima.astro.umd.edu 



CLEAN algorithm and "robust" weighting (robustness 
parameter 0-2) of the visibilities to optimize the signal- 
to-noise and the spatial resolution. Resulting noise levels 
are 0.15 Jy/beam in 0.16 km s^^ channels for the N2H+ 
J = 1 - and HCO+ J = 1 - hnc emission, 0.4, 0.6 
and 1.0 Jy/beam in 0.13 km s^^ channels for the ^^CO 
J = 1 - 0, C^^O ./ 1 - and ^^CO J = 1 - line emis- 
sion respectively. The rms levels in the continuum maps 
are 6.8 and 4.0 mJy beam~^ for the continuum images at 
2.7 and 3mm, respectively. The best naturally-weighted 
resolution is obtained for N2II+ J = 1 — 0, with a beam of 
FWHM of 5.9" X 4.6". Moderate resolution is obtained 
for HCO+ J = 1 - and ^^CO J = 1 - with FWHM of 
7.6" X 5.8" and 9.2" x 5.9" respectively. Due to hmited 
(u, v) coverage, the FWHM of ^^CO J = 1 - and C^^O 
J = 1 — data is 14.3" x 6.7". Integrated-intensity and 
velocity-centroid images were obtained from the cleaned 
spectral-line cubes using a 1 or 2cr clip level. The resolu- 
tion of the continuum images is 9.5" x 6.3" and 6.4" x 
4.9" for 2.7 and 3.3 mm respectively. 

2.2. FCRAO Data 

To obtain information on large spatial scales, we ob- 
served Blc in N2H+ J 1-0 and HCO+ J = 1-0 emis- 
sion with the Five Colleges Radio Astronomical Obser- 
vatory (FCRAO). The data were obtained in very good 
weather. We achieved an rms of 0.08 K (T^) in 24 min- 
utes on source. The N2H+ J = 1-0 and HCO+ J = 1-0 
data were obtained simultaneously using the array SE- 
QUOIA. The beamsize of the FCRAO data is 57.5" at 
93.17378 GHz and 60.5" at 89.188523 GHz. The maps 
cover an area of diameter 12.5' centered on the contin- 
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uum peak of the BIMA array data. 

2.3. Combination of Single Dish and Interferometer 

Data 

For the HCO+ and N2H+ data sets, we were able to 
combine the BIMA array and FCRAO data to obtain 
total power maps of each transition. We uti lized the 
method described bv IStanimirovic et al.l l)1999D to com- 
bine the data in the image plane. 

The FCRAO data were reordered to match the BIMA 
axes (x,y,v) , converted to main beam from antenna tem- 
perature, and rescaled to Janskys. The MIRIAD task 
REGRID was then used to regrid the FCRAO data to 
match the parameters of the BIMA array data cube. 
Slight differences in line frequencies were compensated 
for by a shift in the reference velocity of the FCRAO 
spectra prior to regridding. The shifts required were 
0.885 and 0.016 km s'^ for N2H+ J = 1 - and HCO+ 
J = 1 — spectra, respectively. In the latter case, the 
shift is significantly less than the width of individual 
channels. A map of the FCRAO beam for each transition 
was generated and truncated at the 5% level, creating a 
mask applied to the appropriate single dish data cube. 
The weighting factor for the single dish map was deter- 
mined by a ratio of the beam areas. The composite map 
was then created using a linear combination of the BIMA 
dirty map and the single dish map, followed by decon- 
volution using the combined beam ijStanimirovic et al.l 

3. DUST EMISSION TOWARD BIG 

Figure n shows the continuum detections toward Blc 
at 3.3 and 2.7 mm. The source is located at aj20oo = 
03'^33'^17^878, 5j2ooq = +31°09'31'.'98, which is the peak 
of the 3.3 mm emission. The 2.7 mm emission is poorly 
resolved with positive flux density detected to a radius 
of 13". The 3.3 mm emission is better resolved with the 
core extending to 11". 

The continuum emission is not point-like; Figure ^ 
shows that it is extended to the northwest and the south- 
east of the continuum peak at both wavelengths. The 
continuum emission to the southeast does not coincide 
at the two wavelengths, but in each case it is strong, ex- 
ceeding 4f7 at 3.3 mm at aj20oo = 03''33™18!75, <5j2ooo = 
-F31°Q9'21"1, corresponding to 16.6 mJy beam'^ The 
2.7 mm peak to the southeast also exceeds 3(J (21.7 
mJy beam-i) at a,720oo = 03'^33'"20!20, (5,/2ooo = 
-F31°09'14"5. 

Figure 121 shows the visibility amplitude as a function of 
(it, v) distance at each wavelength. Flux detections above 
the zero-sigma value (dashed lines) are found for only 
a few of the shortest (it, v) distances and are generally 
marginal detections. Based on these plots, we interpret 
the continuum maps as spatially filtered observations of 
a resolved, extended envelope, similarly to the case of 
L483 (see Figure 2 of J0rgensen (2004)). The poorer 
{u,v) sampling at 2.7 mm leads to the recovery of less 
extended emission than we detect at 3.3 mm. 

Using the MIRIAD task IMFIT, we have fit gaussians 
to the 2.7 mm and 3.3 mm maps, fixing the peak flux den- 
sity and position of the peak emission. As Table 13 shows, 
the resultant total flux densities derived are higher than 
those obtained by use of aperture photometry. As ex- 
pected from Figure |21 less flux is recovered at 2.7 mm 



than at 3.3 mm. However, the sizes measured for the 
core are remarkably similar at each wavelength. Based 
on the 3.3 mm deconvolved size, the scale of the detected 
inner envelope is ~ 3100 x 1800 AU at 250 pc. The ef- 
fective radius of this central region (reg = Vab where 
a and b are the major and minor axes, respectively) is 
~ 2400 AU. This is much smaller than the 12000 AU 
effective radius derived from the SCUBA map at 850 /im 
by Kirk. Johnstone & Di Francesco (2006). The mea- 
surement from the JCMT incorporates the entire outer 
envelope of Bl-c, to which our BIMA array data are not 
sensitive. It is normal to measure less flux in interfer- 
ometric maps than in single dish data because interfer- 
ometers preferentially sample structure on small scales. 
Therefore, our continuum estimates only apply to the in- 
ner envelope (i.e., not the region traced by the JCMT). 

At 3.3 mm, it is possible that some continuum emission 
may be due to free- free emission instead of emission from 
dust. We have reduced an archival observation^ from 
the Very Large Array at 1.3 cm in which Blc lies just 
outside the primary beam. There is no detection of 1.3 
cm emission in this map, and the 3cr upper limit is 3.7 
mJy beam""'^. Since the free-free emission is expected to 
be relatively flat, a contribution of this magnitude cannot 
be a signiflcant source of the continuum emission at 3.3 
mm. 

3.1. Mass, Column and Density of the Inner Core 

The mass of the continuum source is easily determined, 
assuming optically thin conditions, from the flux of the 
source and the temperature via the relation 



where is the flux, d is the source distance, is the 
dust opacity and B^(Td) is the Planck function at tem- 
perature Td- Table 13 presents the masses calculated for 
our measured flux densities for different assumptions of 
(3 (and hence k^) and Td- Our continuum maps are spa- 
tially flltered and are not sensitive to the large scale emis- 
sion of the core (i.e., the outer envelope). If this source 
is young, then much of the mass is expected to remain 
in the envelope; therefore, observations sensitive only to 
the inner envelope should be expected to yields masses 
smaller than that measured with a single dish telescope. 

Kirk, Johnstone & Di Francesco (2006) derive a mass 
of 2.4 M0 for Blc based on SCUBA observations at 850 
/im. The clump size (i?oft — 49") and mass were mea- 
sured with CLUMPFIND l)Williams. de Geus fc BliO 
119941) . This estimate of mass and size includes the entire 
outer envelope, and we note that the size scale exceeds 
the maximum scale to which our BIMA array data are 
sensitive: 1.1' at 2.7 mm and 1.3' at 3.3 mm. We con- 
centrate our discussion on the 3.3 mm emission since, as 
discussed above, our flux recovery is better at that wave- 
length. The requirement that the interferometric mass 
estimate be less than the single dish mass estimate al- 
lows us to put some constraints on (3 and Td- The dust 
temperature is > 15 K and [3 is more likely to be closer 

^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc.. 
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Fig. 1. — Continuum emission from Blc at 3.3 mm and 2.7 mm. The respective beams are shown at the bottom right. The 3.3 mm 
contours are plotted at 2, 4, 6 and 8 <t where tr = 4.1 mjy beam"'^. The 2.7 mm contours are at levels of 2, 3, 4, 5, and 6 cr where a = 6.8 
mjy beam~^. 



TABLE 2 
Continuum results 







2.7 mm 


3.3 mm 


Measured Peak Flux Density 


[mjy beam ^1 


47.5 ± 14 


39.0 ± 13 


Total Flux Density (aperture) [mjy] 


71 ± 21 


123 ±37 


Total Flux Density (gaussian 


fit) [mJy] 


89 ± 27 


150 ± 50 


Size ["] 




12.5 X 8.9 


13.8 X 8.8 


PA [°] 




-51 


-30 


Deconvolved size ["] 




10.2 X 1.6 


12.3 X 7.2 


Deconvolved PA [°] 




-66 


-33 



''Measured in apertures of 13" and 11" radius for 2.7 and 3.3 mm, 
respectively. 



40 
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Fig. 2. — Plots of visibility amplitude versus projected baseline 
for our continuum datasets. Points are the amplitudes; the error- 
bars are the formal standard deviations from the mean. The dotted 
line represents the signal expected from the noise. We have only a 
single track each at 109 and 115 GHz, both in the C array. This is 
reflected in the extremely limited (it, v) distance coverage. Those 
bins where the detected amplitude exceeds the value consistent 
with noise are even more limited. The coverage is much better at 
3.3 mm, where we have multiple tracks at 89 and 93 GHz, in both 
the B and C arrays. Similarly, the range of projected baselines in 
which we detect signal well above those consistent with noise are 
also much more extended. 

to 1.0 than to 1.5. For the same dust temperature as 
derived from single dish data (15 K) and (3 = 1.0, the 



mass of the inner envelope is 2.9 ± O.QMq. For Td = 20 
K, the mass is 2.1 ± 0.6Mq. Therefore, the mass of the 
inner core lies in the range 2.1 — 2.9Mq. 

Here, we reiterate that our continuum data are hin- 
dered by excessive spatial filtering. This caveat prevents 
us from rigorously predicting values of /? or dust temper- 
ature in the inner core, since the BIMA flux densities are 
not easily compared either to each other or the existing 
single dish flux density from the JCMT. 

Using the continuum flux, we can also estimate the 
column density of molecular hydrogen within the core 
using the relation: 



N{H2) = 



(2) 



fi rriH Bi,(Td) 

where 5*,^ is the peak flux density, Qm is the main beam 
solid angle in steradians, fi is the mean molecular weight 
(2.33), mn is the mass of atomic hydrogen, k^, is the 
dust opacity per unit mass. Table 13 contains estimates 
of the column density for different values of (3 and Td- At 
3.3 mm, for a temperature of 20 K, the column density 
toward the peak is (6.8 ± 2.0) x lO^^ cm'^, for /3 = 1.0. 
Conversion to extinction using N{}i2)/Av = 10^^ cm^ / 
mag yields a visual extinction of 680 ± 200 mag to the 
central peak of the Blc core. 
Assuming the inner envelope is as deep as it is wide. 
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TABLE 3 

Masses and Column Densities from Continuum Flux Density 



2.7 mm 



3.3 mm 





(5 = 1 


13 = 1.5 


13 = 1 


f3 = 1.5 




0.0049 


0.0034 


0.0039 


0.0025 




1.2 ±0.4 


1.73 ± 0.5 


3.71 ± 1.1 


5.91 ± 1.8 




0.92 ±0.3 


1.32 ± 0.4 


2.86 ±0.9 


4.55 ± 1.4 




0.66 ±0.2 


0.94 ±0.3 


2.06 ±0.6 


3.29 ± 1.0 


cm"-^) 


4.3 ± 1.3 


6.2 ±1.9 


12 ±4 


20 ±6 




3.3 ± 1.1 


4.7± 1.4 


9.5 ±2.9 


15 ± 5 




2.4 ±0.7 


3.4 ±1.0 


6.8 ±2.0 


1.1 ± 0.3 



(cm^ g~i) 
Mass {Td = 12 K) (Mq) 
Mass (Ta = 15 K) (Mq) 
Mass (Ta = 20 K) (Mq) 
Column Density (T^ = 12 K) (lO^^ 
Column Density (T^j = 15 K) (lO^^ 
Column Density (T^j = 20 K) {10^3 



"Derived from Ki, = Kniv/u o)^ where kq is calculated to be 0.01 cm^ g ^ at fo = 231 GHz 
lOssenkopf fc Hennini3ll99l) . 
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Fig. 3.— Maps of integrated intensity in I^CO J = 1 - 0, ^^CO J = 1 - 0, C^^O J = 1 - 0, HCO+ J = 1 - 0, N2H+ J = 1 - 
and SiO J = 2 — 1 toward the source Barnard Ic. (a) N2H+ J = 1 — map shows the double-peaked distribution with peaks offset from 
the 3.3 mm continuum peak of (c), shown as a triangle. The N2H+ data are integrated over 21.2 km s~^ (a range which contains all 7 
hyperfine components); contours range from 2 to 10 it in steps of 2(T where tr = 47 mjy beam~^. Dashed contours show the corresponding 
negative contours, (b) C^^O J = 1 — data clearly trace the core emission and are centrally peaked on the 3.3 mm continuum source peak 
(triangle). The C^^O data are uniform weighted for better resolution (11.7" X 5.4") and integrated over 2.1 km s ^, and the contours 
range from 2 to 4 ct in steps of 1 cr = 0.17 Jy beam^^. (c) SiO J = 2 — 1 data shows no significant detections. Contour levels are 2 and 3cr 
where a = 0.03 Jy beam~^. (d) ^^CO J = 1 — integrated intensity in greyscale and white contours; contours range from 2 to 8 cr where 
(T = 0.16 Jy beam~^. Negative contours of the same range are also plotted. The ^^CO map is integrated over 12.7 km s~^. (e) HCO+ 
map in greyscale and white contours; data are integrated over 17.2 km s~^. The contours range from 2 to 6 cr where cr = 35 mJy beam~^. 
(f) i^CO J = 1 - integrated emission over 4.3 km s ^reveals no emission associated with the continuum peak (triangle). Contours range 
from 2 to 4 cr in steps of 1 cr = 75 mJy beam~^. As in other plots, the dashed contours represent corresponding negative values. 
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then we can crudely estimate the central density by 
adopting the effective diameter as the depth and assum- 
ing a constant density sphere. Then the density is just 
the column density divided by the depth, and the central 
column density is ^ (9.0±2.6) x 10® cm^"^. A comparable 
calculation using 2.7 mm data yields the lower estimate 
of (3.3 ± 1.0) X 10® cm~'^. Reasonable estimates of the 
density therefore range from (3 — 9) x 10® cm~^. 

4. MOLECULAR EMISSION TOWARD BIG 

Maps of the integrated intensity of ^^CO J = 1 — 0, 
i^CO J = 1 - 0, Ci^O J = 1 - 0, HCO+ J = 1 - 0, 
N2H+ J = 1 - 0, SiO J = 2 1 as detected by the BIMA 
array are presented in Figure |3| The outflow from Blc 
is prominently detected in ^^CO J = 1 — and HCO+ 
J = 1 - 0. The N2H+ J = 1 - 0, C^^O J = 1 - clear 
trace the envelope around the central source. The ^^CO 
J = 1 — detections show no emission associated with 
the core, but may trace some compact features in the 
outflow. SiO J = 2 — 1 is not detected. The systematic 
velocity of the source is 6.35 ± 0.02 km s""'^, based on a 
fit to the N2H+ J = 1 — spectrum at the position of 
the continuum source. 

4.1. The Core in N2H+ and C^^ 

N2H+ J = 1 — is strongly detected in Blc, indicating 
dense, cold gas within the core. Figure ^ shows that 
the integrated intensity of N2H+ J = 1 — exhibits a 
double-peaked structure, with both peaks avoiding the 
position of the continuum source (indicated by the tri- 
angle) and the outflow, as shown in Figures |3i andOfc. 
In contrast, Figure Ob shows that the C^^O J = 1 — 
emission is strongly centrally peaked very close to the 
continuum source position. Figure^shows a comparison 
of the distributions of the C^^O and N2H+ emission. The 
emission in a strip along the flattened axis of the core is 
also shown. This figure clearly illustrates the relative po- 
sitions of the peaks in the C^^O J = 1 - and N2H+ 
J = 1 — emission. Despite our significantly poorer 
resolution in C^^O J = 1 — 0, the distribution is strik- 
ingly similar to that observed toward L483 (Figure 10 of 
iJorgenseni (2004)). The relative strength of the emission 
at each N2H+ peak is comparable, which differs from the 
results in L483, in which one peak was noticably brighter 
than the other. 

This morphology is seen across a wide range of young 
pr otostellar objects . It is very similar to that observed 
bv l.TorgensenI l)2004D in the more luminous (9^©) evolved 
Cl ass source L483, and i s also seen in NGC1333 IRAS 
2 ijJorgensen e t al.' zOOla") and sources in the Serpens 
molecular cloud (Hogcrhcijdc 2005). It is also noted in 
the very low luminosity object (VcLLO) IRAM 04191 
(•0.15L^. iLee. Ho fc Whiti^ (2005.) . In Blc, the N2H+ 
emission is not confined merely to the two peaks, but 
extensions are detected to the east and southeast which 
bound the outflow emission in -'^^CO and IICO+. These 
extensions could be the dense edges of a conical cavity 
carved into the envelope by the molecular outflow, as de- 
picted in our schematic diagram of the source (see Figure 
IS)). The projected opening angle of this cavity is ~ 55°, 
based on the lowest contour of the eastern extensions of 
N2H+ emission in the integrated intensity images (Figure 
|3Ji) relative to the position of the continuum source. The 
uncertainty on the opening angle of the cavity is ~ 10°. 



A third extension in N2H+ is observed to the northwest 
along the 3a contours of the IICO+ and CO emission. No 
N2H+ is detected along the other edge of the red-shifted 
outflow emission in the integrated intensity map. 
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Fig. 4. — A comparison between N2H+ and C^**0 emission 
toward Blc. Upper panel: The integrated N2H+ (red contours) 
and C^*0 (greyscale and white contours) emission. The flattened 
direction of the core is nearly north-south and is indicated by the 
dashed Une. The beams are denoted at the lower right. Lower 
panel: The N2H+ (solid line) and C^*0 (dashed line) emission 
along a slice denoted by the dashed line in the upper panel. The 
anti-correlation between the peaks is evident. The ordinate value is 
the ratio of the flux density at the position along the slice relative 
to the peak flux density in the flold. 
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Fig. 5. — A schematic representation of the main features of 
the Blc source. Our interferometric data are insensitive to the 
large scale outer envelope. We detect dense material from the in- 
ner envelope in N2H+ emission, some of which may be conflned 
in a rotating torus centered on the source. The outflow, detected 
in ^^CO and HCO^ has carved cavities in the inner (and presum- 
ably outer) envelope. A centrally heated cavity is indicated by the 
anticorrelation of C^^O and N2H+ at the core center. Dust emis- 
sion remains centrally peaked on the protostar because very high 
temperatures (2000 K) are required for its destruction. 
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Fig. 6. — The moment data for the red and blue CO outflow lobes toward Blc. The greyscale is C^*0 J = 1 — emission from the core; 
the blue contours indicate blue-shifted emission (2.4 - 5.1 km s~^) and the red contours are the red-shifted emission (8.4 - 13.8 km s~^). 
Some red-shifted emission is detected to the southeast of the blue lobe, indicating that the back of the conical cavity is being detected, or 
that the outflow is shifting direction at this location. 

4.2. The Outflow in HCO+ and ^"^CO 

The molecular outflow is clearly seen in both ^^CO 
J = 1 - and HCO+ J = 1 - 0, as shown in Figure EJl 
and We did not detect any SiO J = 2 — 1 emission 
as shown by Figure Et- The detected ^^CO J = 1 - 
and HCO"'" J = 1 — emission exhibits a bipolar mor- 
phology with the blue lobe emission to the southeast 
much stronger than the red lobe emission to the north- 
west. Significantly, some "'^^CO J = 1 — emission is 
also detected coincident with the central core, but HCO+ 
J = 1 — is not detected toward the position of the con- 
tinuum source. As observed in L483 by Jorgensen ct al.'' 
(HomH), the N2H+ data is anticorrelated with the HCO+ 
emission. 

The zeroth moment of ^^CO and C^^O data are com- 
pared in FigureEl The C^^O J = 1 — emission is highly 
centrally peaked on the continuum source and is absent 
along the outflow axis, due to it relatively low abundance. 
The ^^CO J = 1 — outflow appears to have carved out 
the lower density envelope traced byC^^O J ~ 1 — (in- 
tegrated over the central velocities associated with the 
source) in a similar manner as the very dense gas traced 
by N2H+ J = 1 — in Figure Ot- In fact, the optically 
thin species C^^O J = 1 — traces well both the contin- 
uum peak and small scale structure within the inner en- 
velope, detected on the edge of the ^^CO J = 1 — emis- 
sion along both the blue and red lobe s . Thi s is consistent 
with the findings of lArce fc Sargent (2006) in their sur- 
vey of outflow sources with the Caltech Millimeter Array. 
The orientation of the outflow from Blc is approximately 
—55° east of north. Based on the blue lobe, the opening 
angle appears to be ~ 35° in projection. Figure also 
clearly indicates the presence of red-shifted CO emission 
at the leading edge of the blue lobe. 

4.2.1. Spitzer Mid-infrared Emission 

Figure [7| shows the outflow emission from the source 
Blc as imaged by Spitzer with IRAC at 4.5 /im 



(jJorgensen et al.ll2006j) and the ^^CO outflow emission 
from the BIMA array. The IRAC 4.5 /im band has been 
found to be a strong tracer of outflows which is likely 
due to the presence of H2 pure rotational transitions and 
the C O fundamental vibrational m ode within the band- 
pass l)Noriega-CresDO et al.ll200l) . The outflow clearly 
extends for several arcminutes on either side of the driv- 
ing source, Blc. The 4.5 /xm emission is a strong tracer 
of molecular hydrogen, typically associated with proto- 
stellar jets. An "S-shaped" morphology is evident in the 
Spitzer map. Such S-shaped jets are interpreted as pre- 
cession of the jet because of the p resence of an (unseen) 
binary (e.g., iHodapp et al.ll2005|) . Near the source, the 
outflow appears quite symmetric; however, at larger dis- 
tances, the blue lobe widens while the red lobe appears 
to either split into two separate sequences of "bullets" 
or be confused with an outflow from a different source. 
The projected linear extent of the outflow detected by 
Spitzer is dependent on which distance one adopts to 
the Barnard 1 cloud. At 250 pc, the 6' extent of the 
blue lobe from the central source indicates a distance of 
- 90000 AU, or 0.44 pc. 

Comparison of the Spitzer data with the BIMA ar- 
ray data illustrates that the 4.5 /im emission lies along 
the central axis of the CO outflow very near the driving 
source. Like the molecular hydrogen emission, there is a 
bend in the ^^CO emission, located precisely where the 
CO emission transitions from blue-shifted to red-shifted 
emission. We are constrained by the single pointing of 
the BIMA observations to a single primary beam of cov- 
erage in the CO emission. More extended ^^CO J = 1 — 
measurements are needed to compare the morphology of 
the molecular hydrogen (the driving jet) to the morphol- 
ogy of the entrained gas mapped by the ^^CO. It is in- 
teresting that the position of the Spitzer emission peaks 
are anti-correlated with the CO peaks in the BIMA data, 
which is likely due to the difference in excitation condi- 
tions between the H2 transitions that sample warm (few 
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hundred K) gas along the jet shocks and ^^CO J = 1 — 
which preferentially traces entrained, cold (~ 10 K) gas. 

4.3. 13 CO J = 1 - Emission 

Figure |2t shows the integrated emission detected from 
I'^CO J = 1 — 0. In contrast to the other CO isotopo- 
logues, the ^^CO data do not clearly trace the core or the 
outflow. This is likely due to missing short spacings in 
the interferometric map which limit the sensitivity of this 
observation both to structure and emission. The peaks 
detected in ^^CO may be associated with the outflow or 
the walls of a cavity carved by the outflow. 

4.4. Moment Maps of the Core and Outflow 

We produced moment maps of the combined N2H+ 
J = 1 — emission from the FCRAO and BIMA array 
datasets over a velocity range of 1.57 km s"^, or 10 chan- 
nels. After creating the moment maps, we masked out all 
data values in the moment zero and moment one maps 
at positions with values less than 0.24 Jy beam"^ km 
s~i in the moment zero map, which is approximately 1.5 
times the rms level per channel, 0.14 Jy beam^^. 

Figure |H1 shows the first and zeroth moments of the 
combined N2H+ J = 1 — emission. The velocity gra- 
dient suggests that the molecular core is rotating about 
an axis aligned with the outflow. As in the integrated 
intensity map from the BIMA array of Figure 13^, the 
zeroth moment map of the combined dataset illustrates 
a double-peaked distribution with peaks offset from the 
position of the continuum emission. This morphology 
is expected from the projection of a torus of dense gas 
surrounding the core center (see Figure 13)). We similarly 
interpret this double-peaked feature as a rotating torus 
with the blue-shifted emission lying predominantly to the 
north and the red-shifted emission lying to the south as 
indicated by fits to the spectra and the features of the 
first moment map. Figure |H1 shows the moment one map 
in greyscale and the moment zero map in contours. We 
do not have sufficient resolution in our N2H"'' J = 1 — 
spectra to separate out all the hyperfine splitting (hfs) 
components in the spectra. The presence of the hfs struc- 
ture complicates the interpretation of the velocity field 
from the line emission, which can be alleviated by using 
the isolated component. The moment map is taken over 
the isolated line component of N2H"'" J = 1 — Fi F 
01 - F12 at a v elocity offset by - 8.0064 km s^^ from 
the central line l|Dore et al.ll2003j) . The velocity range 
over which the moments are taken is 1.57 km s""'^ (10 
channels). 

The moment map of N2H+ shows the same extensions 
to the east and northwest as the integrated intensity 
map shown in Figure Ot- Both the eastern extensions 
detected in N2H+ are red-shifted compared to the adja- 
cent gas, with the south eastern one containing the most 
red-shifted gas in the core. The northwestern extension 
appears blue-shifted. A faint extension is detected to the 
west of the core which could be the southern edge of the 
cavity created by the red-shifted lobe of the outflow. 

For comparison, in Figure we show the zeroth mo- 
ment map of the FCRAO data alone over the whole 
mapped area. The Blc source is part of a broader 
distribution of dense N2H"'" J = 1 — emission which 
forms a bridge in emission to the s ources to the south : 
Bla (IRAS 0330H-3057) and Bib ijHirano et al.llT997|) . 



The distribution is singly peaked, not surprisingly at 
the position of the brightest peak detected in the BIMA 
and combined maps, as shown in Figure El through im- 
ages of the BIMA array data and combined data con- 
volved to the FCRAO beam size of 57.5". We note 
that the N2H''" emission is comparable in extent to the 
850 jiim SCUBA continuum scale of the c ore (diameter 
1.6' lKirk. Johnstone fc Di Francescoll2006D . In contrast, 
BIMA emission alone is tracing scales < 1.3'. 

The axis of rotation of the combined map of Figure |H1 
lies at ~ —45° east of north, which is close to the —55° 
orientation angle of the outflow (see S I4.2|I . indicated by 
the dot-dashed line. However, the axis of rotation is not 
orthogonal to the plane of the N2II+ peaks, which lies at 
^ 10° and is indicated in Figure|Slby a dashed line. The 
offset between the plane of a potential torus (as estimated 
from the double peaks) and the rotation axis as derived 
from the flrst moment map is 55°, not 90°, as might be 
expected. Similarly, the offset between the torus and the 
outflow axis is apparently 65°. 

The absence of N2II+ emission along the outflow axis 
is likely due to the carving out of a cavity by the jet as- 
sociated with the outflow. The combined FCRAO and 
BIMA zeroth moment maps of N2H+ and IICO+, shown 
in Figure 1101 mirror the morphology observed in the 
BIMA data alone. The outflow emission traced by HCO"*" 
J = 1 — neatly fills the cavity in the N2II+ emission. 
The most strongly blue shifted emission of IICO+ does 
not lie at the largest distances from the central source. 
Similar to FigurelHl there is evidence for redshifted emis- 
sion at the edge of the blue lobe; the outer edge is cer- 
tainly less blue-shifted than the material closer to Blc. 

4.5. Kinematics of the Core and Outflow 
4.5.1. Position- Velocity Diagrams 

We have taken a position-velocity cut along the plane 
of the torus (through the continuum peak position and 
the two N2H+ J = 1 — peaks) as sh own in Figur e 
|S| using the Karma program "kpvslice" ('Gooch"1996"). 
Figure ITTl shows the PV emission of the isolated N2II+ 
,7=1 — component. The emission is continuous across 
the molecular core, but emission peaks are well separated 
in velocity space. They are symmetric about the velocity 
of the source for this HFS component (—1.64 km s"^). 
The distribution of the emission is very narrow in velocity 
space, limited to only a few velocity channels. Nonethe- 
less, there is a clear velocity gradient across the core with 
a discernible shift in the centroid velocity across the two 
off-center peaks. 

Figure El shows a position- velocity slice taken along 
the outflow axis through the continuum peak. The orien- 
tation of the slice is shown as a dashed line on FigurelTUI 
The distribution of the outflow in velocity space is rela- 
tively confined in the HCO"^ J = 1 — emission (which 
contains both single dish and interferometric data). The 
velocity features are more extensive in the ^^CO J = 1 — 
emission, even though it is limited to the spatial frequen- 
cies detectable to the interferometer. The red-shifted 
material seen at the leading edge of the blue lobe is ev- 
ident in the PV diagram and no blue-shifted emission 
is observed at offsets beyond that boundary. The red- 
shifted material is detected only in the ^^CO J = 1 — 
emission. It is not seen in HCO+ ,7=1-0, although 
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Fig. 7. — IRAC Spitzer data at 4.5 /xm from |j0rgenseii et al.l 120061) are compared to moment maps of the ^^CO J = 1 — emission over 
the blue and red lobes. The correlation between the position angle and central axis of the outflow is excellent between the ncar-IR and 
millimeter data. The inset shows the larger scale extent of the molecular hydrogen emission seen from Spitzer. ^■^CO J = 1 ~ contours 



are from 3-10 Jy beam ^ km s ^, in steps of 1 Jy beam 
circle indicates primary beam of a BIMA antenna. 
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Fig. 8. — First (greyscale) and zeroth (contour) moment maps 
from the combined FCRAO and BIMA array observations of N2H+ 
1-0 emission. The contours are intervals of 0.5o" from 1.5 to 6 
a whore a = 0.115 Jy beam~^ km s~^. The cross marks the 
position of the continuum peak. These maps are taken over the 
single isolated component of the N2H+ hyperfine transitions, Fi F 
= 0,1 — > 1,2. The white dashed line represents th e sl ice taken to 
produce the position vs. velocity diagram of Figure lTTl and defines 
the orientation of the "torus" at 10° east of north. The dot-dashed 
line represents the position angle of the outflow (-55° east of north) 
as me asur ed from IRAC 4.5 fira and BIMA ^^CO J = 1 - data 
(see 

Figure ITUl shows that the strongest blue shifted emission 
is detected close to Blc. At larger distances from the 
source, the HCO"'" J = 1 — emission is dominated by 
velocities closer to that of the source. 

4.5.2. Spectra 

Figure [T^ shows the distribution of spectra (spatially 
binned to an area comparable to the beam size) across 



The cross marks the continuum peak of the BIMA data. The large 

the core in the combined FCRAO and BIMA data. Fits 
to the hyperfine components of the N2H"'" spectra were 
done using the CLA SS software's "hfs" fitting routine 
ijBuisson et al.' '2002) with up-to-date weights and fre- 
quencies Pore et al. 200^. These fits were used to 
determine LSR velocities (Vlsr), intrinsic line widths 
(Aw), total optical depths (T(toi)), and excitation tem- 
peratures (Tox). The distribution of line strengths is very 
atypical for N2H^. The hyperfine ratio suggests that the 
central velocity components should be strongest when 
compared to higher and lower clusters of hyperfine com- 
ponents. In Blc, the central component is often the same 
strength as the higher velocity component, which is in- 
dicative of high opacity or self-absorption of the brightest 
component. Solutions to the hfs fitting were not sub- 
stantially improved by using only the outer triplet and 
the isolated component of the seven hyperfine lines to 
produce fits to the data for the LSR velocity and the 
linewidth; therefore, the fits shown in Figure ll^SI are all 
for the fits derived with all seven components. 

All seven hyperfine components from the combined 
data were fit simultaneously in CLASS. We spatially 
binned the data onto a grid with spacing 5" in R.A., 
and 6" in declination. The fits reflect centroid velocities 
ranging from 5.94 to 6.68 km s~^. The range of realistic 
fine widths (FWHM) is ~ 0.3 to 1.5 km 8"^ The source 
velocities are systematically bluer to the north and red- 
der to the south as has been discussed in § 14.41 

FigurelT^shows the distribution of Vlsr, AF and line 
opacity (sum of the peak optical depths for all seven hy- 
perfine components) across the core on the same grid 
as shown in Figure El The gradient in Vlsr. is well 
defined across the field. Higher values of AV are ob- 
served along the outflow axis where the cavity has been 
carved in N2H+ emission than across the central core 
where linewidths are quite uniform. Some of the best 
fit solutions (assuming a single excitation temperature 
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Fig. 9. — Comparison of the zeroth moment N2H+ J = 1 — emission on large scales from FCRAO, the BIMA array and combined 
FCRAO and BIMA array data, where the latter two maps have been convolved to the FCRAO beam size of 57.5", shown at the top left 
of the FCRAO image. These maps are taken over all seven hyperfine components. The contours of the FCRAO and combined maps are 
at intervals of IOct from 20 to 100 cr where a = 1.7 Jy beam~^ km s~^. The BIMA array map has substantially less flux than either the 
FCRAO or combined maps. The contours plotted range from 4 to 12 Jy beam~^ km s~^, in steps of 4 Jy beam~^ km s~^. A cross marks 
the position of the 3 mm source. The BIMA array data, convolved to a similar beam, peaks at the same position as the FCRAO data 
alone. 



S 31 09 00 




Fig. 10. — Upper panel: Moment zero maps of HCO+ (greyscale) 
and N2H+ (contours) show the displacement between the HCO+ 
and N2H+ emission. The HCO+ emission is confined to a cavity 
carved into the dense gas by the outflow. Contour levels are 2 to 6 
cr in steps of 0.5 cr where cr = 0.11 Jy beam~^ km . The dashed 
line shows the position of th e position-velocity slice taken through 
the outflow (see Figure [T^ . Lower Panel: The first (greyscale) 
and zeroth (black contours)moment maps of HCO"*" 1-0 emission 
from the Blc outflow. In the combined map of FCRAO and BIMA 
data, the blue lobe of the outflow remains brighter than the red 
lobe. The white contours show the moment zero N2H+ emission. 

for all components) shown in Figure 1131 produce very 
large total optical depths for the northern and south- 
ern emission peaks. These values indicate that at least 
one component is very optically thick and all compo- 
nents may be optically thick. Attempts to fit the spectra 
with lower, fixed estimates of r(tot) do an increasingly 
poor job of fitting the heights of the components. The 
software's minimum measurable optical depth is 0.1; in 
these instances, the emission must be optically thin, but 
we cannot discriminate between values of T[tot) lower 
than 0.1. The maximum value yielded is 30, obtained 
when components have equal strength. 

Figure E| shows only the isolated component of the 
N2H+ spectra, with a velocity range limited to -4 to 
km s""'^. These spectra reveal differences both in the peak 
of the lines and in line shape across the core, with some 
positions exhibiting a dip in the center of the line profile. 
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Fig. 11. — Position- velocity diagram of N2H"'' emission along the 
plane of the dual peaks observed in the moment map of Figure |S] 
Only the isolated component of the hyperfine transitions is used. 
The central position is coincident with the continuum peak. The 
systemic velocity of the core for the HFS component is also marked 
with a dashed line. The velocity gradient across the core is evident. 



TABLE 4 

Parameters Derived from N2H+ hfs fitting 





Mean 


rms 


Min 


Max 


Visr [km s 1] 


6.31 


2.34 


5.94 


6.75 


AV [km s-l] 


0.61 


0.37 


0.26 


1.85 


T(tot) 


6.99 


8.37 


0.100 


30.0 



^The maximum and minimum solutions for T(tot) 
based on the hfs routine are 0.1 and 30 respectively. 
Values of 0.1 indicate all components are optically 
thin, while 30 indicates that the component strengths 
are equal. 

which could indicate that even this line may be optically 
thick. It is also obvious that regions of poorer S/N pro- 
duce poorer fits, higher (likely unphysical) estimates of 



5. DISCUSSION 
5.1. Column Density and Chemistry 
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Fig. 12. — The position-velocity diagram of HCO+ J = 1 - 
and ^-^CO J = 1 — emission produced along the outflow axis. The 
greyscale shows the HCO+ J = 1 — emission; '^^CO J = 1 — 
emission (BIMA only) is shown in contours of 20 to 90 percentile 
in steps of 10 percent. The slice is centered on the position of the 
Blc co ntin uum source. The orientation of the slice is indicated on 
Figure 021 

5.1.1. Column Density of N2H^ 

In order to estimate the N2H"'" column density, we 
require estimates of the FWHM of the hnes and esti- 
mates of the Hne opacity. Figure El shows the distribu- 
tion of hne widths and opacities derived from the N2H+ 
J — I ~ spectra. We derived the N2H+ column den- 
sity using a simple curve of growth determined by as- 
suming LTE excitation conditions at = 12 K, a line 
width of 0.7 5 km s~^, and the RADEX escape proba- 
bility code^ l|Schoier &: van dei~Ta k ,2005) . This simple 
approach takes into account line opacity, but a more de- 
tailed envelope model with proper density and tempera- 
ture slope will be required to get more accurate estimates 
(Matthews et al. in preparation). Figure [TBI shows the 
column density distribution of N2H+ J = 1 — emission 
in Blc. Since the LTE analysis predicts opacities on the 
order of unity while the spectra and hfs component fit- 
ting suggest significantly higher opacities, it is likely that 
the excitation is not in LTE (i.e., Tex < 12 K). However, 
Figure lTCl is illustrative as it shows the same depression in 
N2H^ emission seen in the zeroth moment map of Figure 

[13 

5.1.2. Column Density of C^^O 

To derive the density at the core center, we estimate 
the column density at the peak of C'^^O J = 1 — emis- 
sion. We can derive the column density from the expres- 
sion 



N, 



3.34 X 10^^ 



ly ^'^ (1 — exp(— /iz^/fcTcx)) 
/ TuBdV T 



exp(— Jft,i^/2fcrcx) 1 ~ exp(— r) 



(3) 



where v is the frequency in GHz (109.78217 GHz), fi is 
the dipole moment in Debye (0.11 Db), Tex is the exci- 
tation temperature, and J is the lower rotational level 
for the transition. We assume Tkin = Tex which is con- 
stant for all rotational levels, LTE is valid for the C^'^O 

3 www.strw.leidenuniv.nl/ ~moldata 



gas and that the optical depth effects are completely ac- 
counted for by the factor r/(l — e^'^). We evaluate the 
integral over the line width from the peak of the moment 
zero map, integrated over a linewidth of 1 km s~^. We 
estimate r from the peak temperature of the line via the 
expression 



hv 



r cak « — 1 - CXp(-T) 

^'^'^^ k exp(ft,z//fcrkin — 1) 
where Tkin is the kinetic temperature. From the spec- 
trum at the peak position, we evaluate r for Tpcak = 2 
K and T^n ~ TLx = 12 K. We find t = 0.244, which is 
consistent with the C^^O emission being optically thin. 
The strength of the emission at the peak is 1.16 Jy 
beam"""^ km s""'^. Using the conversion factor for our 
BIMA array synthesized beam at 109 GHz (1.12 K/Jy 
beam~^), J TuBdV is 1.3 K km s~^. Substitution of 
these values in equation (j^J yields a column density of 
1.03 X 10^^ cm^^ (the uncertainty, based largely on the 
flux calibration, is ^ 30%). Utilizing the abundance ra- 
tio [C ^^01/[H2l 1.7x 10-^ (Frerking. Langcr & Wilso^ 
Il982t) yields an H2 column density of 1.75 x 10^^ cm~^. 
The effective radius of the beam in the C^^O observa- 
tion is 1.8 X 10^® cm. Assuming the core is as deep as 
it is wide, the density at the peak is ^ 9.5 x 10^ cm~'^. 
This value is three times the estimated central density of 
3. 2 X 10^ cm~^ from the analysis of dust emission maps 
bv lKirk. Johnstone fc Di Francescol l)200fiD , but a magni- 
tude less than the estimate based on 3.3 mm continuum 
data in S I3.1I The difference is partly attributable to the 
different spatial filtering between the observations. The 
2.7 mm data, with more comparable {u,v) coverage to 
the C^^O data, yields densities of (3.3± 1.0) x 10^ cm"^, 
closer to the value derived above. 

Aside from the uncertainties introduced by the differ- 
ent degrees of spatial filtering, the discrepency between 
the H2 column density measured with dust and C^^O 
most likely arises due to depletion from standard ISM 
abundances, leading to underestimates of iV(H2). Al- 
though the grains must be heated in Blc due to the 
centrally peaked distribution of C^^O, some of the gas 
could remain frozen onto the dust grains. Assuming 
the dust represents the total column of the inner enve- 
lope and the gas is depleted, then the [C^®0]/[H2] ra- 
tio is not 1.7 X 10"^ but instead 1.3 x 10"^ (from 2.7 
mm dust continuum), a factor of 13 less. This value 
is similar to the values measured in Class sources by 
iJorgensen. Schoier &: van Dishoecklllmil . 



5.2. A Depression in N2H^ Emission near the Core 

Center 

As demonstrated by the integrated intensity and mo- 
ment maps, the N2H"'' emission is diminished near the 
core's center. We are able to produce a more coher- 
ent picture of the central "cavity" in N2H"'' emission 
by reweighting the N2H+ data presented in Figure |S1 to 
uniform when creating the maps. The downweighting 
of shorter baselines degrades the S/N of the individual 
channels significantly, but the moment map produced re- 
tains high signal by utilizing all seven hyperfine compo- 
nents. Figure ITTl shows the moment zero map produced 
from the uniformly weighted BIMA array data combined 
with FCRAO data. A high threshold (3cr) was placed 
on inclusion of individual channels into the moment map 
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Fig. 13. — N2H+ spectra across the Blc inner core. The spectra are integrated across the core in grids of 5" X 6", which is slightly 
greater than the beam area. Therefore, the spectra are independently sampling the kinematics of the core. The velocity range is -10 to 20 
km s~^, and the intensity range is -0.5 to 2 Jy beam~^. Contours of the moment map are plotted below the spectra. Contours range 
from 0.35 to 0.91 Jy beam~^ km s~^ in steps of 0.14 Jy beam~^ km s~^. The red lines show the fits to the spectra from CLASS and the 
fitted values for Vlsr E'nd a-rs shown at the left and right corners, respectively, of each grid square. 



to avoid channels between components. The cross marks 
the position of the continuum source, and the triangle 
represents th e pointing c enter of the JCMT observations 
presented in IMatthews fc Wilson ( 2002). The ring de- 
picts the 14" beam of the JCMT. The scale of the evap- 
orated cavity is such that it would not have been resolved 
by the JCMT. The scale of the cavity is 2.7" x 2.15", 
at a position angle of 7° . This corresponds to 675 x 540 
AU at a distance of 250 pc. 

The depression of N2H"'" emission in the inner envelope 
and along the outflow path is likely due to the destruc- 
tion of N2H+ by molecules that are released from grains, 
most notably CO and H2O. Chemical models of an evolv- 
ing core with a central luminosity source (i.e., Class 0) 
suggest that by 10^ years a cavity of ~ 700 AU will be 
seen in the N2H+ abundance and emission ijLee et al.l 
|2004|) . This is the direct result of stellar heating providing 
enough energy to evaporate CO from grains within the 
core center. This model may not be directly applicable to 
Blc due to differences in masses and/or luminosity, but 
it provides a rough estimate of the time that has passed 
since stellar ignition. In addition, observations and mod- 
els suggest t hat NgH+ will also be destroyed w ithin the 
outflow ijBachiller et alJ2001l:rBergin et al.ll998(l . and we 



see a strong anticorrelation between the position of the 
outflow and the presence of N2H+ emission. 

A central depression in N2H+ emission has also been 
detected in the transition Class O/I source L483 by 
Udrgensenl 1200^ and in the l o w lum in osity object IRAM 
04191 by IBelloche fc ATdTl ll200l . ILee. Ho fc Whit^ 
()2005(1 estimate the inner edge of the IRAM 04191 cavity 
has a mean radius of 1400 AU. This is a factor of four 
times larger than the cavity detected in Blc. In a sim- 
ple model where cavities grow over time, this suggests 
that Blc could be even younger than IRAM 04191's es- 
timated age of 3 X 10"^ yr since the onset of accretion 
|Andre, Motte & BacmannI Il999fl . The gradients and 
line widths measured in N2H+ J = 1 — are strikingly 
similar in these two cores, indicating strong dynamical 
similarities, although we stress that IRAM 04191 is one 
of the new Class of very low luminosity objects (VeL- 
LOs, Di Francesco et al. 2006; Young et al. 2004) while 
Blc is a relatively luminous (4 ± 2 L©) Class object 
(Matthews et al. 2006, in prep). 

We can estimate a kinematic age for this source based 
on the outflow data. The projected linear extent of the 
emission coupled with the distance from the source yields 
estimates of the source age. Using the brightest peak of 
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Fig. 14. — Distribution of derived values of VlS-Ri linewidth and total optical depth based on the fits shown in Figure ITsl The gradient 
across the source is obvious in the derived values of V^sR- The linewidths are broader along the direction of the outflow and are narrower 
along the plane of the torus detected in integrated intensity. The optical depths show a large range of values from optically thin to very 
optically thick {>> 1), indica ting that several hfs components may be optically thick. Table 3 shows the mean values of the parameters 
derived from the fits to Figure [T31 

the blue lobe detected in ^^CO J = 1 — and assuming an 
inclination angle of 45°, the dynamical age is 3.7 x IC* yr, 
comparable to that of the source IRAM 04191 in Taurus 
(3 X 10'* years). To produce an age on the order of 10^ 
years, the outflow inclination angle must be ~ 15° (i.e., 
the outflow must lie almost along the line of sight). 

The absence of N2H+ in IRAM 04191 is at- 
tributed to depletion within the cold core interior 
(jBelloche. Andr6. Despois fc Blind er 2002). The strong 
centrally peaked detection of C^^O J = 1 — emission 
from Blc (see Figure 0) shows that depletion is not sig- 
nificant in Blc, and the anticorrelation between N2H+ 
with the core center and outflow lobes suggests that 
N2H+ is absent due to destruction by CO. In addition, 
the strong outflow detected with Spitzer (Figure [Tj) is 
highly collimated, suggestive of a young source, but its 
spatial extent suggests that the collapse of the driving 
source may not have been as recent as that in IRAM 
04191. These characteristics are similar to those of L483 
which is thought to be in transition to a Class I object 
lITafalla et al.ll200ni) . 

ILee. Ho fc Whitd l|2005D note that the centers of 
the molecular outflow and envelope surrounding IRAM 
04191 are offset from the the continuum source by ap- 
proximately 560 AU. They suggest that the two cavities 
seen to the south of that source could be due to a binary 
in which one companion is too young to have any appre- 
ciable dust emission. There is also evidence for binarity 
within Blc, particularly in the behavior of its outflow 
(S I4.2|l . which shows evidence of precession. The cavity 
center is also offset from the continuum source position. 
There is as yet no detection of a second dust peak nor a 
second outflow in Blc. 



5.3. Rotational Support 

The observed global rotation velocities within the Bl 
cloud are insufflcient to su pport the cloud again st col- 
lapse by a factor of ~ 8 l|Bachiller et al.llT990j) . The 
ages of embedded but optically visible ob jects LkHa 327 
and L kHg 328 are between 4-6 xlOfyr jCohen fc Kuhl 
119791) . Based on this, Bac hiller et al.l l)1990D conclude 
that another mechanism must be providing substantial 
support to the Bl cloud. Within Blc, the two peaks of 
the rotating torus can be used to derive the rotational 
energy relative to the Vlsr at the core center of Blc 
itself. 

The inner edges of the torus are offset from the posi- 
tion of the continuum peak. For a sphere, the rotational 
energy is given by 



Er, 



pVvt 



(5) 



where p is the mean density, V is the volume, and Vj. is 
the radial velocity relative to the center. Based on es- 
timates of the central density with C-'^^O (^ | 5. 1.211 and 
dust emission (jKirk. Johnstone fc Di Francesco .2006^) , 
we adopt a mean density of 5 x 10^ cm~^ for Blc. 

The radial velocity of the two peaks relative each other 
is ^ 0.06 km s~* (as derived from Figure [TS|l . Thus, we 
estimate the radial velocity relative to the core center 
to be a relatively small shift of ^ 0.03 km s~*. The 
mean offset of the peaks from the core's velocity center 
is ~ 4 ", or 1000 AU (for d — 250 pc), encompassing a 
volume of 1.4 x 10*^ cm'^. Thus, based on equation (j^J, 
t he rot ational energy o f the c ore is ^ 1.2 x lO'^^ ergs. 

iMatthews fc WilsonI l)2002() derived the energetics of 
the Blc core based on estimates of the magnetic field 
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Fig. 15. — N2H+ spectra of the isolated hyperfine component of N2H+ across the Blc core. The spectra are integrated across the core 
in grids of 5" X 6", which is slightly greater than BW sampling. The velocity range is -4 to km s~^, and the intensity range is -0.25 
to 1 Jy beam~^. Contours of the moment map are plotted below the spectra. Contours range from 0.35 to 0.91 Jy beam~^ km s ~-^ in 
steps of 0.14 Jy beam~^ km s~^. The fits shown as red lines are derived using CLASS for all seven hyperfine components fFigurc ll3i . The 
fitted values for Vi^sR S'^'d AV are shown at the left and right corners, respectively, of each grid square. 



dispersion and found rough equipartition between the 
gravitational, kinetic and magnetic energies (all ^ 10^* 
ergs). The rotational energy is thus likely negligible to 
the support of the core compared to the other means of 
support available through nonthermal gas motions and 
magnetic field support. 

The radius of the torus 1000 AU) detected in 
Blc leads naturally to a discussion of a "pseudo-disk" 
(j^^alli fc Shu 1993) around the protostar within this core. 
Such a disk is magnetically suppo rted and grows with 
time. iMatthews fc Wilson! l)2002t) estimated the mag- 
netic field strength in the low density regime around the 
Bl cores to be on the order of 30 fiG. While the field 
strength within the Blc core may be higher, there is no 
Zeeman splitting measurement of sufficient spatial reso- 
lution to isolate the field strength of the individual cores 
within Bl. Assuming the prese nce of a pseudo-disk as 
described bv lCalli fc ShTl we can use their equa- 

tion (24) to estimate the age of the disk around Blc. For 
a field strength of 30 /iG and a sound speed of 0.35 km 
s~^, a disk of radius 1000 AU has an age of 1.8 x 10^ yr. 



consistent with the expected age of a Class protostar. 
6. SUMMARY 

Utilizing interferometric and single dish millimeter 
data from the BIMA array and FCRAO, we have probed 
the structure of the envelope of the source Barnard Ic 
in Perseus. We detect a powerful molecular outflow in 
12 CO J = 1 - and HCO+ J = 1 - which has carved a 
conical cavity into the envelope of Blc, destroying N2H+. 
The prominence of the outflow driven by this source is 
also evident in Spitzer IRAC data recently published by 
iJorgensen et all l(2006|) . Comparison of the CO molec- 
ular outflow and 4.5 fim data reveals that positions of 
changes in jet direction are mirrored by changes in ve- 
locity in the CO emission. 

We have detected significant N2H+ J = 1 — emission 
from Blc. Data from FCRAO reveal N2H+ emission is 
of comparable scale to single dish dust continuum maps. 
High resolution data (created by uniform weighting of 
the BIMA data combined with the FCRAO data) reveal 
clearly a small cavity in which dense gas has been de- 
stroyed, most likely by heating by the outfiowing gas. 
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Fig. 16. — Column density distribution of N2H+ J = 1 — 
emission in Blc. Contours are sum of the N2H"'" J = 1 — emission 
over all channels with signal exceeding 3 times the rms level. 
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Fig. 17. — Moment zero map of the N2H"'" emission from Blc 
with uniform weighting. The new resolution is 3.09 " X 2.52", 
almost a factor of 2 improved over the data in § 14.41 The cavity 
can now be seen as an absence of N2H+ emission offset to the west 
of the continuum peak, marked by a cross. The triangle marks the 
position of the source as determined from the SCUBA data; the 
ring d epict s the JCMT beam of 14". Contours are the same as in 
Figure [TSl 

We have detected a gradient in Vlsb. of N2H+ consistent 
with rotation; moment maps reveal evidence of a rotat- 
ing torus offset from the continuum peak position. Blc 
is similar to the more luminous source L483. Both L483 
and Blc show C^^O toward the source peak and an anti- 
correlation of C^^O with N2H"'", indicating that the inte- 
rior is warm and that CO has been desorbed from grains 
in that core interiors. In a subsequent paper (Matthews 
et al. 2006, in preparation), we will model the centrally 
heated region within the Blc core using new data from 
the SubMillimeter Array. 

Based on the 3.3 mm dust continuum emission, we 
estimate that the mass of the Blc inner envelope is in 
the range of 2.1 — 2.9Af0. This is consistent with the 
mass of 2.4 Mq derived from the JCMT observations 



ijKirk. Johnstone &: Di Francescdl2006fl for a mean dust 
temperature of 15 K. However, the 3.3 mm measurement 
is only an estimate of the inner envelope due to the spa- 
tial filtering of the BIMA array data. 

It is interesting to note that this core exhibits many 
features of the standard picture for an evolving young 
stellar object. In continuum dust emission, the core has 
a dense centrally peaked interior. However, the N2H+ 
J = 1— emission appears elongated along an orientation 
roughly orthogonal to its outflow. The two N2H+ peaks 
suggest the presence of a rotating torus (with r ~ 1000 
AU) with a central cavity carved out by the outflow and 
heating by the central source. We argue the heating has 
released ^^CO and its isotopes from grains near the cen- 
tral source, destroying N2H+. N2H+ survives in the rest 
of the inner envelope, as indicated by the torus (or pseu- 
dodisk) and the N2H"'" emission which brackets the out- 
flow as traced in ^^CO J = 1 - and HCO+ J = 1 - 0. 
As is the case in L483, the double peaks of N2H+ (torus) 
do not lie orthogonal to the outflow orientation. In Blc, 
the offset is ~ 65°. 

Existing polarimetry data l)Matthews fc Wilsonl l2002). 
interpreted in the case of a constant field orientation 
through the core would indicate a magnetic field direc- 
tion roughly parallel to the outfiow orientation. As yet, 
there is no estimate of the field strength within Blc it- 
self; howe ver, using the estima te from the surrounding 
Bl cloud llMatthews fc Wilsonll200^ and the model of 
iCalli ShiJ l)1993t) . the source ag-e is estimated to have 
a lower limit of 1.8 x 10^ yr. The dynamical age estimate 
of 3.7 X 10* yr is highly dependent on the inclination of 
the outflow. 

In our forthcoming paper (Matthews et al., in prepa- 
ration), we will present higher transition observations of 
^^CO and C^^O from the Blc outfiow, as well as measure- 
ments of the 1.3 mm continuum from the SubMillimeter 
Array. With existing data, and new data from IRS on 
Spitzer, we will generate the spectral energy distribution 
and use it to constrain models of T(r) and p{r) which was 
not possible based on the weak continuum detections re- 
ported here. 
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